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Abstract. The Kovacs or crossover effect is one of the peculiar behaviours exhibited 
by glasses and other complex, slowly relaxing systems. Roughly it consists in the non- 
monotonic relaxation to its equilibrium value of a macroscopic property of a system 
evolving at constant temperature, when starting from a non-equilibrium state. Here, 
this effect is investigated for general systems whose dynamics is described by a master 
equation. To carry out a detailed analysis, the limit of small perturbations in which 
linear response theory applies is considered. It is shown that, under very general 
conditions, the observed experimental features of the Kovacs effect are recovered. The 
results are particularized for a very simple model, a two-level system with dynamical 
disorder. An explicit analytical expression for its non-monotonic relaxation function is 
obtained, showing a resonant-like behaviour when the dependence on the temperature 
is investigated. 
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1. Introduction 

Glassy behaviour has been quite an active field of investigation in recent years. A 
review of the main, experimentally observed features of real glasses, as well as of several 
models showing dynamical properties analogous to them can be found in (TJ E]. In 
relaxation experiments, linear response functions usually exhibit non-exponential decay. 
In particular, a Kohlrausch- Williams- Watts (KWW) law is often found. Moreover, a 
laboratory glass transition, in which the properties characterizing the macroscopic state 
of the system become frozen, is also observed. The transition is associated to a fast 
increase of the relaxation times upon decreasing the temperature. Along reheating, 
hysteresis effects show up, with the system going back to equilibrium following a curve 
which differs from the equilibrium one. In these reheating experiments, the curve 
overshoots equilibrium, and the difference between the actual value of the property 
of interest and its equilibrium value as a function of the temperature shows a non- 
monotonic behaviour. This is a typical example of a memory effect: the behaviour of 
the system depends on its entire thermal history, and not only on the instantaneous, 
initial value of the property under study. 

One of the first and simplest experiments revealing a memory effect was designed 
by Kovacs [U [5]. Polyvinyl acetate was equilibrated at a high temperature T , and 
the "direct", monotonic, relaxation function (p(t) of the volume to a low temperature 
T < T was measured. Then, after equilibrating the system again at T , it was rapidly 
quenched, this time to a lower temperature T\ < T, at which it isothermally relaxed for 
a time t w . This time t w was not enough for the system to reach its equilibrium state at 
temperature Ti. Then, the temperature was abruptly increased to T. The time t w was 
chosen such that the volume at t w was equal to the equilibrium value at T. The observed 
behaviour of the system for t > t w turned out to be quite peculiar. The volume did not 
remain constant, but it first increased, passed through a maximum at a certain time 
t k > t w , relaxing afterwards to the equilibrium value at T. The results of the Kovacs 
experiment are qualitatively sketched in figure HJ although the property considered there 
is the energy E, instead of the volume V, because the former is the quantity that will 
be studied in this paper. The non-monotonic behaviour for t > t w is represented by 
the curve K(t), while 4>E{t) describes the direct relaxation of the energy from To to T. 
In the Kovacs experiment, the pressure P was kept fixed along all the processes. Then 
the observed behaviour means that the knowledge of the state variables, P, V, and T, 
does not fully characterize the macroscopic state of the system, since the subsequent 
evolution is different for systems with the same values of the state variables but with 
different thermal history. This is the reason why this experiment is said to show a 
memory effect. Besides, this non-monotonic behaviour, sometimes called the Kovacs 
"hump" , displays some characteristic features. For fixed initial and final temperatures 
T and T, the magnitude of the maximum increases as the intermediate temperature Ti 
decreases, i.e. the maximum is higher the larger the second temperature jump T — T\. 
Besides, the position of the maximum moves to the left, in the sense that t k — t w is 



The Kovacs effect: a master equation analysis 



<E> 




<E>(V 



AT 



<e^(T) 



AT' 



<E%(TJ 



Figure 1. Schematic representation of the Kovacs experiment described in the text. 
The dashed green curve 4>e (t) represents the direct relaxation from To to T. The dotted 
red curve stands for the part of the relaxation from To to Ti, which is interrupted by 
the second temperature jump, changing abruptly the temperature from T\ to T at 
t = t w . After this second jump, the system follows the solid blue curve K(t), which 
reaches a maximum for t = tk and, afterwards, approaches 4>E{t) for very long times. 



a decreasing function of T — Ti. Finally, for very long times, when the Kovacs hump 
function K{t) decreases towards the equilibrium value at T, it approaches the direct 
relaxation curve <f>E- 

The Kovacs effect has been investigated in many systems showing slow glassy 
dynamics. Kovacs himself presented a phenomenological description in [5]. Later on, it 
has been analyzed both analytically and numerically in several models. This includes 
a molecular dynamics analysis of a realistic model of ortho-terphenyl [6], simulations 
of a 3d spin-glass |7J, of classical and quantum p-spin models [5], and of kinetically 
constrained models (9j \TU\ . Also, analytical investigations have been carried out in the 
context of simple models. Examples are the Glauber Id Ising model [11], the 2d XY 
model [12], domain growth and trap models [IB], the parking lot model for granular 
compaction [TJ], the HOSS models for a fragile glass [15], and a distribution of two-level 
systems [16]. All the above systems exhibit a behaviour resembling the experimental 
Kovacs "hump" , with a maximum whose height increases with the temperature jump 
T ' — Ti. The time tk for the maximum usually decreases with the jump T — Ti, although 
in some models tk — t w appears to be almost independent of the jump and, also, the 
tendency towards the direct quench curve is not clear [TU [16] . 

In this work, we would like to improve the current understanding of the Kovacs 
effect by trying to identify the origin of its generality. In particular, one of the main 
aims of this work is to prove that it can be considered as a general feature of systems 
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whose dynamics is described by a master equation, i.e. homogenous Markov systems. 
Some of the models mentioned above belong to this kind of systems [101 [TU [151 EE] ■ 
Besides, we will assume that the transition rates verify detailed balance, although the 
results can be also valid in some cases in which this condition is not fulfilled, as discussed 
in the last section of the paper. In order to carry out a detailed analytical analysis, the 
temperature jumps will be considered small enough, so that linear response theory can 
be applied. It will be shown that, in this limit, the Kovacs effect as well as its main 
features are still present, clearly indicating that it is not a non-linear effect. 

The plan of the paper is as follows. Section [2] is devoted to a general discussion of 
the linear response to a temperature jump of a system whose dynamics is described by a 
master equation. The obtained results are applied to the Kovacs experiment in section [3j 
and a general expression for the Kovacs "hump" function K(t) is derived. Some general 
features of this function are discussed in section HI In section [5J the two level system 
with dynamical disorder is introduced, and some of its properties are shortly reviewed. 
The Kovacs effect in this system is investigated in section |6j The simplicity of the model 
allows us to obtain an explicit analytic expression for K(t), which is analyzed in the 
light of the general results obtained previously. Section [7] contains a summary of some 
of the main results and conclusions of the present work. Finally, some calculations are 
presented in the two appendices. 

2. Linear response to a small temperature jump 

Consider a system whose dynamics is described by a homogenous Markov process. The 
conditional probability Tt-f(i\j) of finding the system in state % at time t, given that it 
was in state j at time t' < t, obeys the master equation 



Here Wj& is the (time-independent) transition rate from state k to state i. The above 
equation has to be solved with the initial condition T (i\j) = 5ij. The probability Vi(t) 
that the system be in state i at time t is obtained from the conditional probability and 
some initial condition 7^(0) through 
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The set of both Tt-t'(i\j) and Vi{t) completely determines the Markov process. From 
(JTJ) and fl2]) it follows that Vi(t) also obeys the master equation, 
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It is convenient to introduce the matrix W defined by [19] 
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so that ([3]) can be written in matrix notation as 

±V{t) = WP(t) , (5) 

where V(t) is the column matrix with elements Vi(t). The formal solution of this 
equation is 

p(t) = e m V(0) (6) 
or, in component notation, 

w = E( eCT %Pi(°)- w 

3 

Equivalently, the solution of (GQ) is 

T t - t ,(i\j) = (e^) ■ (8) 



The transition rates often depend on a set of macroscopic parameters X = 
{Xi, X 2 , . . . , X r }, such as the temperature, pressure, density, external fields, and so 
on. In some cases, the values of the above parameters can be externally controlled. 
Along this paper, attention will be restricted to this sort of situations. Suppose that 
a system is evolving being X^ ^ = x!f \ . . . ,Xr°^ the values of the macroscopic 

parameters and = W(X^) the corresponding evolution matrix. The probability 
distribution V^°'(t) will obey the master equation, 

=W (o)^(o) (t)- (9) 

Next, consider the evolution of the same system, but with parameters given by 
X = X(°) + AX. The new transition matrix can be decomposed as 

W = W(X) = W (0) +W, (10) 

where, 

W = W(X (0) + AX) - W(X (0) ) . (11) 

The new probability distribution V(t) will obey (jSJ), and it can be decomposed in the 
form 

Vit) =V {0 \t)+V'(t). (12) 

Obviously, both W and V vanish for AX = 0. Substitution of (HQJ) and (JT2J into © 
and use of yields 

^-V'(t) = W^V'(t) + W'P (0) (t) + W'V'Ct). (13) 
at 

In the linear response approximation, the limit AX — ?- is considered, and it is assumed 
that the last term on the right hand side of (TLB"]) can be neglected, since it is at least of 
second order in the deviations AX. 

In the following, the above results will be particularized for an idealized kind of 
experiments. Starting from a given initial condition V(0), the system evolves for a time 
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t w with values X^°) of the macroscopic parameters. Then, at t = t w the parameters are 
instantaneously changed to X = X^ ) + AX, being AX very small. The probability 
distribution of the system for t > t w is given by ( Tl2l) and ( 1131) . Integration of the latter 
after neglecting the nonlinear term leads to 

V\t) = f dt'e^ wi0) W'P^(t') , (14) 



valid for t >t w , while V'(t) = for t < t w . In component notation, 

Jt w j k 

t >t w , with 

= (e^ w<0) ) . (16) 
The average value of property A at time t > t w , will be denoted by 
{A{t\t w )) = Y, A ^)^ 

i 

where A{ is the value of the property A in state i. From (fl2|) and (fl5j) it follows that, 
again for t >t w , 

(A(t\t w )) = (A(t))<® + /" d*'itU(M') , (18) 



with 

(A(t))(°) = ^^ 0) (t) (19) 

i 

and 

R A(ft') = ^^^^rS,(.|j)w;. fc ^ 0) (t'). (20) 



Equation ([18]) implies that 

A [(A(t\t w )} - (A(t))<®] = -R A (ft w ). (21) 

The above relations indicate that R,A{t,t w ) can be interpreted as the linear response 
function in the property A of the system to a modification of the external parameters 
X. It is worth stressing that the system has not been assumed to be at equilibrium at 
time t = t w , when the perturbation is switched on. Suppose now that this is the case, 
and that the equilibrium distribution of the system has the canonical form 

Ve,i(T) = -^=r, (22) 

where /3 = (ksT)" 1 , ks being the Boltzmann constant and T the temperature. 
Moreover, Ei is the energy of the system in state % and Z(T) the partition function 

Z(T) = J2 e ~^ 1 - ( 23 ) 
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In addition, let us suppose that the only external parameter that is changed at 
t = t w is the temperature, which is modified from to + AT or, equivalently, 
/?(°) = (k B T^Y l is changed into 0® + A/3 with A/3 = -(^T^j^AT. Under these 
circumstances, the response function is given by 

R e>A (t, = - = A/3^(A(t)£(t'))i 0) , (24) 

t > t' > t w . Here (A(t)B(t')}^ denotes the equilibrium correlation function at 
temperature T^°\ 

i 3 



A proof of f}24"]) is sketched in Appendix A 



3. The Kovacs experiment 

Let us analyze the experiment proposed by Kovacs jH |5| and qualitatively depicted in 
figure [TJ To characterize the non-monotonic behaviour of the average energy, let us 
decompose it for t > t w in the form 

(E(t\t w )) = (E) e (T) + K{t), (26) 

where the Kovacs "hump" K(t) satisfies 

K{t) > Vt , lim K[t) = lim K(t) = . (27) 

t—>t w t— »oo 

The specific shape of the hump depends both on the final and intermediate temperatures 
of the system, T and T 1; respectively. It exhibits a maximum at a certain time For 
fixed T, the height of the maximum increases with the temperature jump T—Ti, while its 
position relative to t w , % — t w , decreases. For long enough times, (E(t\t w )) approaches 
the curve 4>E{t), corresponding to the relaxation at temperature T starting from the 
equilibrium state at temperature T . 

To render the problem tractable by analytical methods, the temperature jumps 
will be assumed to be small, i.e. |AT| 1 and |AT'| 1, where AT = T — To and 
AT' = T\ — T. Then, the results derived in the previous section can be used. For the 
time interval < t < t w , the system relaxes at temperature T\ = T + AT' starting from 
an equilibrium situation at T^°\ and direct application of (118j) and (124"1) gives 

<JE?(f)> = (E) e (T - AT) - (A/3 + A/3') jf dt' ^{E(t)E(t')) e , (28) 
< t < t w , with 

AT AT' 

In fl28l) . the equilibrium correlation function (E(t)E(t')) e must be evaluated at 
temperature T, since calculating it at temperature T would introduce nonlinear 
corrections in the temperature jumps, which are neglected in the linear response theory 
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we are using. In Appendix B it is shown that, for systems in which detailed balance 
holds, it is 

(E(t)E(t')) e = {E)l + ^(A^) 2 e-^ W , (30) 

a 

where a is the index for the eigenvectors of the matrix W, — \( a > is the eigenvalue 
corresponding to a, and 

AE^ = Y,EM a) - (31) 

i 

Here ip^ is the i-component of the eigenvector a. The prime to the right of the 
summation symbol in (130]) indicates that the equilibrium eigenvalue is excluded. Thus, 
the equilibrium time autocorrelation function of the energy is a monotonic decreasing 
function of the time interval t—f, since it is a linear combination of decaying exponentials 
with all its coefficients being positive. Particularization of ( l30l) for t = t' provides an 
interpretation for AE^ a \ 

al^(E") e -(E)l = Y^(AE^r, (32) 

a 

i.e. (AE^) 2 is the contribution of the a-th mode to the equilibrium energy dispersion. 
Equation fl30|) can be expressed as 

(E(t)E(t')) e -(E)l = aU E (t-t'). (33) 



with 



and 



>> E (t) = J^b^e- X(a)t (34) 



&<«> - {AEia))2 (35) 
The weights have the properties 

> 0, Y^ bH = 1 • ( 36 ) 

a 

This relaxation function is "normalized", in the sense that 4>e{0) = 1, and it 

verifies that 0£(oo) = for all temperatures. Returning to (T2"8"|) . and using f l3"2"j) and 



(E(t)) = (E) e (T) + o%Ap - (A/3 + A/3 r y E [1 - <p E (t)} 

= (E) e (T + AT') + (A/3 + A0')(T%<f> E (t) ■ (37) 

This expression shows the monotonic relaxation of the average energy towards its 
equilibrium value at temperature Ti = T + AT'. Now, the condition fixing the time t w 
at which the temperature is increased is employed, 

(E(t w )) = (E) e (T) (38) 
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i.e., using ( 157)) . 

Mt w ) = a n t~~ a = , 'r:,' , = * • < 39 i 



A/3' | 


AT' 




A/3 + A/3' 


AT 


1 + 1 


AT'| 



The last equality defines the variable x taking values in the interval < x < 1. 
Accordingly with the Kovacs experiment described above, at t — t w the temperature is 
instantaneously raised from Ti to T = T\ — AT'. To write the expression of the average 
energy for t > t„.( 1T8l) will be used again, but in this case the reference evolution is 
provided by (157)) . which describes the energy relaxation at temperature T = T + AT'. 
Consequently, the response function R E {t,t') to be used corresponds, in principle, to 
a non-equilibrium situation. However, since Re is proportional to W, it follows that, 
in the linear approximation, the distribution function in ( 120)) can be replaced by the 
equilibrium distribution at temperature T. That means that the equilibrium response 
function has to be considered again. In this way, it is easily obtained that 

(E(t\t w )} = (E} e (T + AT') + (A/3 + A/3')4<M*) 

+ A/3' f &t>^(E{t)E{t')) e , (40) 

and, taking into account again (155)1 . 

(E(t\t w )} = (E) e (T) + (A/3 + A/3')a 2 E <f) E (t) ~ A/?Vl^(t - t w ) . (41) 

This equation is similar to the one obtained by Kovacs by means of a phenomenological 
theory [5]. The difference is that, in the phenomenological theory, the arguments of the 
response functions are reduced times, in the spirit of Narayanaswami-Moynihan-Tool 
theory of supercooled liquids [ZD [22j [25J , instead of actual times as in ()4T)) . 

To characterize the evolution of the system for t > t w , i.e. after the second 
temperature jump, it is convenient to use a dimensionless Kovacs hump (compare with 
<m>) by 

(E(t\t w )) - (E) e (T) ^ (E(t\t w )) - (E) e (T) 
{ ] cr|A/3 " (E) e (T - AT) - (E) e (T) ' 1 ] 

This is the function which is analyzed in Kovacs' experiments [U [5]. Although it is 
not explicitly shown in our notation, K*(t) is actually a function of both t and the 
"waiting" time t w . Therefore, K*(t w , t) will also be used for the Kovacs hump when the 
dependence on t w should need to be emphasized. Substitution of ( HTI) into ( 142)) gives 
K *u\ = <pE(t)-x<p E (t-t w ) ^ 
1 — x 

where x = <t>E(f w ) is the function of the temperature jumps defined by (1591) . 

In the next section, some general properties of this function will be derived. They 
follow from the form of the normalized energy autocorrelation function 4>E{t), as given 
by (154]) . For the analysis, it is useful to introduce the functions 

D,( ()s -*M.W = -g|, (44) 

D2(t) =-* tH - Mt)]= -m. (45) 
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In the above expressions the prime is used to indicate derivative with respect to time. 
Both Di(t) and D 2 (t) are strictly positive decreasing functions of t, as a consequence of 
4>E{t) being a linear combination of exponentials with positive coefficients. This implies 
that D\(t) must tend to a well defined limit for t — > oo, 

lim D x (t) = Df* > . (46) 

t— >oo 

If 7^ 0, the relaxation is basically exponential for very long times. This is 

compatible with exhibiting non-exponential behaviour in most of its relaxation, 

since the values of (f> E (t) for which the exponential decay is observed can be very 
small, say < 0.01. This is the case, for instance, of the energy relaxation at 

low temperatures in the one dimensional Ising model with Glauber dynamics [T7[ IT8] . 
On the other hand, if = 0, the relaxation is non-exponential even for very large 
times. The above comments follow directly from the definition of D\{t) in f )4*4|) . 

4. Properties of the relative response function K*(t) 
4-1. It is bounded between and <pE{t) for all t >t w 

To begin with, it will be shown that for those times t\ such that K*(t\) = 0, it is 
K*'(ti) > 0, i.e. K*(t) never crosses the time axis with a negative slope. From ( j43l) . 



1 — X 

If K*{t\) = 0, equation ( )43i) implies that x = 4>E(ti) / (f>E{ti — t w ), and substitution of 
this into (I4T1) leads to 

K *,r h) = Mtl)[Dl(tl-t w )- gl(ti)] _ (4g) 

1 — X 

Since the function Di(t) is a monotonically decreasing function of t, the term inside the 
square brackets is positive and, therefore, 

K*'(ti)>0, (49) 

as indicated above. Moreover, the equality sign only holds if 4>e{Px) — 0, i.e. for t\ — > oo. 
As K*(t) vanishes by definition at t — t w , it follows that K*(t) > for all t > t w , tending 
to zero in the limit t — > oo. In addition, K*(t) has an upper bound that follows directly 
from (1451) . by taking into account that < </>e(£ — t w ), 

K\t)<(j> E {t). (50) 

In the experiments [HE], it has been observed that K*(t) < <pE(t — t w ). As 
(^Eif) < (pE(t — t w ), we have derived here a more restrictive inequality, perhaps as a 
consequence of having restricted ourselves to the linear response regime. 
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4-2. There is only one maximum of K*(t) 

As K*(t) is a regular positive function of t, vanishing both at t = t w and for t — > oo, 
it must exhibit at least a maximum for t > t w . It will be proved that actually there is 
only one. Suppose that a stationary point occurs at t = %, so that 

1 — x 

where (l4"7j) has been used. Equation (lo"Tj) defines as a function of i.e. of the ratio 
of the temperature jumps |AT|/|AT"|. The time tk will correspond to a maximum or a 
minimum of K*(t) depending on the sign of the second derivative at tk, 

K *"^ = 73^ " ^Efa - • ( 52 ) 

(|5Tj) gives that x = 4> E {tk) / 4>' E {t k — t w ), which leads to 

K*"(t k ) = [D2 (t k _ ^) _ £, 3 ( tjfc )] j (53) 

where we have introduced the function I?2(^) defined in fH5l) . The term in brackets is 
again strictly positive, since D 2 (t) is also a positive monotonically decreasing function 
of t. Therefore, taking into account that 4>e is a monotonic decreasing function of time, 
it follows that 

K*"(t k )<0, (54) 

and tk must correspond to a maximum. Moreover K*"(t k ) = is only possible if 
4>' E (tk) = 0) i- e - f° r tk OO) when K*(t) approaches zero with horizontal tangent. It is 
also clear that there can only be one maximum of K* (t) since, on account of continuity, 
between two maximums there should be at least a minimum. 



4.3. Behaviour of K*(t) for \AT'\ > |AT| 

In addition to the limits |AT| 1 and |AT"| 1, required by the linear analysis 
developed above, here the case in which |AT| |AX"| will be addressed. In this limit, 
equation ([55]) gives 

4>E(tw) = x ~ 1 - e, (55) 

where e = |AT|/|AT'|. This means that 0^(t) is still very close to unity for t = t w , i.e. 
t w is much smaller than the characteristic relaxation time of the energy at temperature 
T. The Kovacs hump, given by (j4"51) . can be written as a function of s — t — t w , t w , and 

e, as 

K%s)=(l + 1 -) Ms + tw )-^± 



e 



+ 705s(«) + O (t w ) + O (f ) • (56) 
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As a function of s, K* has a has a maximum at s = Sk that, keeping only the first two 
terms on the right hand side of (156)) . is given by the solution of the equation 



D 2 (s k ) = — ^ • (57) 

6 7/ 1 67, 



1 u; 



The above expression implies that Sk is an increasing function of t w or, equivalently, it is 
a decreasing function of |AT'|, for AT fixed. This can be seen by taking the derivative 
with respect to t w in both sides of (1571) to get 

w/ ^gfc _ 0gfa) - 1 - t w (j)' E (t w ) 

The right hand side of this equation is negative, because the numerator is a strictly 
decreasing function of t w vanishing for t w — > 0. Moreover, D 2 was shown to be a 
monotonically decreasing function of time, following that dsk/dt w > 0. 

Let us analyze equation ( l56l) deeper. Two rather different scenarios can arise, 
depending on the behaviour of t w with e, when the latter is very small. Consider first 
that to lowest order it is t w ~ 0(e). For very short times, it is 

<l> E (t) - 1 - £>i(0)*> ( 59 ) 
where it has been used that -Di(O) = — 0'(O). Particularization of the above relation for 
t — t w and use of fl55l) yields 

*■ = D^Y (60) 
Substitution of this value into f )56|) leads to 

K%s) = Ms) + j$fi +0(e). (61) 

Therefore, when t w ~ 0(e), the maximum of the Kovacs function is located at a value 
sjk verifying 

^ 2 (4 0) ) = f ^A(O). (62) 

Consequently, tends to a finite, well defined value for t w — > 0. This is consistent with 
the behaviour observed in some models, where s^ seems to be almost independent of 
the magnitude of the second temperature jump [TU ITS] . 

Consider next the case t w e. This may happen in systems where there is a time 
window for which t 3> [Di(0)] _1 , so that fl59|) does not apply, but is still very close 

to unity and ( 1551) can be accomplished inside that window. A well known simple model 
exhibiting this behaviour is the one-dimensional Ising model with Glauber dynamics [20] 
in the low temperature region [T71 [18] . A typical time dependence in this time window 
is 

<p E {t) ~ 1 - At\ (63) 

where A > is a constant and 7 a real parameter in the interval < 7 < 1. Then, 
equation ( 1551) gives 

«« = Q) V7 (64) 
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and equation (I56"j) takes the form 

K*(s) = Ms) + 0(e- 1+ ^y (65) 

The above discussion supports the following scenario: as |AT'| increases, K*(s) 
approaches the relaxation curve 4>e{s) because the maximum of the former Sk moves to 
smaller times. It is interesting to note that the stretched exponential 

<pKww(t) = e~^\ (66) 

with r being a characteristic relaxation time and < 7 < 1, which is often used to fit 
experimental data of the relaxation of supercooled liquids and other complex systems 
[HE], has the short time behaviour given by (163"]) . 

44. Behaviour of K*(t) for \AT'\ < |AT| 

For I AT' I < I AT|, (EHD leads to 

I AT' I x , , 

m=— x « 1 - < 67) 

i.e. x = 4>E(t w ) <C 1 and, therefore, t w 3> 1. As the equilibrium energies at the 
temperatures T and T + AT' = T\ are much closer than those for T and T — AT = To, 
the curve describing the relaxation of the energy at temperature T\ crosses (E) e (T) at 
a later stage of the relaxation, as compared with the case |AT| <C |AT'| . 

The function Di(t), introduced in (144p . tends to a well defined limit D^, as 
indicated in fj4*6|) . Therefore, D[(t) vanishes as t — > 00 and this allows the approximation 
<f> E (tw + s) = e ^<t>E(U)-D 1 (t w )s-±D> 1 (t w y+- ^ xe -Di(t w )s ^ ^ 6g ^ 

expected to be valid for D[(t w )s 2 <C 1. Thus, equation (T4"3]) gives 

K*(s) ~ — *— \x e - Dl(t ™ )s - x <Ms)l - x \ e - Dl ^ )s - (f> E (s)] . (69) 
1 — x 

This function exhibits a maximum at s = Sk verifying 

4>' E {s k ) = -D x {t w )e- D ^ s * . (70) 

Two possibilities must be considered at this point, corresponding to a vanishing and 
no n- vanishing long time limit, of Di(t), respectively (see (I46p ). If = 0, it is 

D\{t w ) <C 1 and equation ( ITDjl can be approximated by 

<t>' E {s k )~- Dl {t w ) = ^^ (71) 

as long as Di(t w )sk <C 1, something to be checked a posteriori. Since 4> E (t w ) <C 1, it is 
concluded from the above relation that 4>' E (sk) 3> 4>' E (t w ) and hence s k <C t w . Therefore, 
Sk increases as t w increases and it diverges for t w — > 00, but remaining always much 
smaller than t w , i.e. t w s k 1. Moreover, equation ( J7T1) also implies that 

SkD^ty,) ~ -Sfe0E(sfc), (72) 



The Kovacs effect: a master equation analysis 



14 



and the right hand side of this relation is very small since 4>e{s) tends to zero when s 
goes to infinity. This proves the consistency of the assumed dominant balance used to 
solve ( 1701) . The maximum value of the Kovacs function is obtained by substituting Sk 
into pi}, 

K* max = K\s k ) = x [ e ~ D ^ - Ms k )] ~ x . (73) 

Summarizing: for |AX"| <C |AT| and = 0, the position of the maximum Sk 

decreases as |AT'| increases, while the height of the maximum ~ 4>E(t w ) has the 

opposite behaviour. This agrees with the experimental observations by Kovacs [HE]. 
It is worth stressing that the stretched exponential ( 1SS1) . verifies that D\{t) — » for 
t — > oo, and therefore it fits into the case just discussed. 

Consider next that > 0, so that equation (!69l) reduces to 



K*(s) = x 



-D 



e 1 -<Pe(s) 



■ xip(s) , (74) 

where (f(s) = exp[— £)^°°^s] — 0e(s). Therefore, i^*(s) factorizes into a function of t w 
times a function of s. This means that the maximum of K*(s) will occur at a time 
such that ip'(s) vanishes, and it is independent of t w , 



(°o)> 



-D^e-^^ , (75) 



i.e. s k °°^ is the formal limit of s k given by ( 170|) for t w — > oo. The height of the maximum 
in this case is 



x 



max 



(76) 



It is again a decreasing function of the waiting time t w , or an increasing function of 
the second temperature jump |AT'|, in agreement with the experimental observations. 
On the other hand, the position of the maximum, s^, now does not diverge but tends 
to a finite limit. Also this is consistent with the behaviour found in some models, in 
which the position of the maximum of the hump seems to be independent of the second 
temperature jump [TH US] . 

4-5. Asymptotic behaviour of K*(t) for long times 
Let us start our analysis from (|4"3I) written in the form 

K*(s) = Ms) + lMt w + s)- Ms)} ■ (77) 
For long enough times, 4>E(t w + s) can be approximated by 

Mt w + s) = Ms)e~ Dl{s)tw , (78) 

valid for D[(s)t 2 w < 1. Remember that D[(s) -> for s -> oo. Use of fl7S]) into (1771) 
gives 

K(s) ~ Ms) + ^ [e~ D ^ - 1] , (79) 
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for s ^> 1. As in the previous section, two cases must be analyzed separately. If 
-Di(s) in the long time limit, the above expression simplifies to 

K(s) ~ Ms) , (80) 

i.e., the Kovacs function K*(s) approaches the relaxation function </>e(s) for long times, 
such that Di(s)t w <C 1. This is analogous to the experimental observation jU [5]. On 
the other hand, if the logarithmic derivative does not vanish in the long time limit but 
D x {t) d[°° } > 0, (GSD takes the form 

K*(s)~Ms) 1 T ( :f w) (81) 

1 - <PE[tw) 

or 

In K(s) ~ In Ms) + hi - - |^ , (82) 

for s ^> 1. Consequently, for long enough times plots of In K(s) versus s corresponding to 
different values of t w can be collapsed on ln0^(s), by subtracting an adequate constant 
quantity for each value of t w . 

5. A two-level system with dynamical disorder 

Here the general scenario developed in the previous sections will be particularized for a 
simple model, perhaps the simplest one exhibiting the Kovacs effect. It is a two- level 
system (TLS) with dynamical disorder. There are two possible states of the system 
that will be denoted by 1 and 2, respectively. The difference of energy between the 
states is e, and the energy barrier between them, measured from the excited state 1, is 
V. This barrier V is not fixed but fluctuates in time between two values, V + and V^, 
being V + > V-. These fluctuations are described by a dichotomic Markov process with 
constant transition rate 7. A sketch of the model is presented in figure [3 
The transition rates between states 1 and 2 are given by 

W21 = W(l ->• 2) ± = ae~ pv± = v± , (83a) 
Wf 2 = W{2 -> 1)* = ae- p{v±+e) = v±e~ pe , (836) 

where a is a characteristic attempt rate to cross the barrier and the ± signs correspond 
to the two possible values of the energy barrier. Let V u (i,t) denote the probability of 
finding the system in state i — 1, 2 at time t, being in addition the height of the barrier 
V u ; v = ±. The master equation for the model reads 

lp(t) = WV(t) (84) 



with 



V 



r+(2,t) 

P-(M) 



55) 
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and 



/ -0+ + 7) 



W 



\ 



7 




7 

-{v + e-P E + 7) 

-(y- + i) 

7 v_ 





7 



-/3e 



-Pe 



(86) 



+ 7) / 



The equilibrium distribution V e is the stationary solution of (J8 



V, 



2(1 + P) 



1 

V 1 / 



(87) 



where we have introduced the parameter 



Note that the matrix W in ( tool) satisfies the detailed balance condition. The eigenvalues 
— A^ and eigenvectors ipw of W are the solutions of the equation 

W^ (a) = -A (a V (a) • (89) 

In addition to the null eigenvalue corresponding to the equilibrium distribution, there 
are other three eigenvalues given by 

A (1) = 7 + z/ M (l + p)-K, (90) 

A (2) =-f + u M (l + p) + K, (91) 

A (3) = 2 7 . (92) 
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The associated eigenvectors are 



P 



7 



2(1 + p) 2 v /[ K + 5(l + p)] 2 + 7 2 



(2) 



7 



(3) 



2(1 + P? v /[ K -5(l + p )]2 + 7 2 



1 -1 

2(1 + P) P 

V 1 / 



/ K + S(l+p) \ 

/ 7 \ 

K+(5(l+p) 

7 
-1 



V i 

/ K-6(l+p) \ 

/ 7 I 

<5(l+p)-K 

7 
-1 

V i / 



(93) 



(94) 



(95) 



In the above expressions, the positive parameters 
V- + v+ 



V M = 

6 = 



2 

V- — Vj, 



K= vV + ^l + p)' 



(96) 

(97) 
(98) 



have been introduced. 

For any two matrices V\ and P 2 defining states of the system, a scalar product is 
defined as [T5J 

i 



i=l 



7*1,^2, 
Pe.i 



(99) 



Then, the eigenvectors ip^ verify the orthogonality condition 

(v> (Q) ,V> (a,) ) =<W- (100) 

and they constitute an orthonormal basis for the solutions of the master equation, 



v(t) = v e + J2 a{a) ^ {a)e ~ x(a)t 



(101) 



a=l 



The coefficients are given by 

a (a) = (^W ; p(o)) . (102) 

In this way, the dynamics of the system is completely solved, since the evolution of the 
probability distribution V(t) has been obtained for an arbitrary initial condition V(0). 

The marginal probability V u (t) of finding the TLS with a barrier V u , v = ±, at time 
t, regardless of the state of the system, is given by 

P„(t) = P„(l,t)+P v (2,t). (103) 
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Similarly, 

V(r,t) = V+(r,t) +V-(r,t) (104) 

is the probability of finding the system in state r = 1,2, without taking into account 
the value of the energy barrier. It is interesting to realize that 

a (3) = (^(3)^(0)) = - V+ (l, 0) - V+{2, 0) + 0) + V-{2, 0) 

= P_(0)-7>+(0), (105) 

i.e. a® = if the initial probability is symmetrically distributed between the two 
possible values of the barrier. Thus the third eigenvalue A'- 3 -' is associated with the 
relaxation towards configurations in which both values of V have the same probability. 

An expression for the equilibrium correlation function of the energy can be written 
down by particularizing fl30|) . 

(E(t)E(t')) e = (E)l + (A£^)V A(Q) ^') , (106) 

a=l 

with AE^ given by ( 13~TT) . The term corresponding to a = 3 has been omitted since it 
is AE^ = 0. Taking the origin of energies at the ground state 2, it is 

AE^ = J2 e i4 a) = e \4 a) + 4 a) ] (107) 



2 



;i08a) 



or, explicitly, 

AE (D = _ £ / P k + *(1 + p)+7 

V2(l + p) 2 v /[ K + 5(i + p )]2 + 7 

^ = - e,[~ r ~ Z^±B±M2 , (1086) 
V2(l + p) 2 ^/[ K _5(i + p )]2 + 7 2' 

For the equilibrium energy dispersion it is found 

4=(A^)%(A^)) 2 = II ^. (109) 



The relaxation function of the energy follows by applying fl34|) . 

Mt)=b^e~ xWt + b^e- x{2)t , (110) 

with 

2[ K + <5(l + p)] 2 + 7 2 ' V ; 

l[- w + tf(l + p) + 7 ]» 

2 [k-5(1 + p)] 2 + 7 2 1 ; 

The long time behaviour of 0e(£) is given by 

B (t) ~& (1) e" A(1 \ t^oo, (112) 
since A^ < A*- 2 * 1 . Therefore, the long time limit of the function Di(t) defined in (jUJ) is 
= A (1) > . (113) 
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This is a relevant feature of the present model since, as discussed in section |3J the 
long time limit of D\(f) controls some of the key properties of the function K*(t w ,s) 
characterizing the Kovacs effect. 

In the rapidly fluctuating barrier limit defined by 7 ^> 5(1 + p), the system is 
equivalent to a TLS with an effective barrier, 

A (1) ~ u M (l + p), A (2) ~2 7 + z/ A f(l + p), (114a) 
6 (1) ~1, 6 (2) < 1. (1146) 

Therefore, the relaxation approaches an exponential decay as 7 increases. 

In the slowly fluctuating barrier limit, 7 <C 5(1 + p), the static disorder case is 
recovered. The results tend to those of an ensemble of two kinds of TLS having the 
activation energy e but different barrier heights, V + and VL, respectively, 

A (1) ~ u+(l + p), (115a) 
A (2) ~ + (1156) 



Besides, 



6« ~ - , 6 (2) ~ - . (116) 



6. The Kovacs experiment in the disordered TLS 

Next, the results of submitting a disordered TLS to the Kovacs experiment will be 
discussed. In the linear response approximation, valid for small temperature jumps, the 
function K*(t) given in (T43l) characterizes the evolution of the energy after the quench 
from the intermediate temperature T\ to the final temperature T, at time t w . The latter 
is determined by ( 1391) . that using (lllOp becomes 

b m e -^t w + b (2) e -xvh w = l&T'l s x ; (117) 



which gives t w as a function of the temperature jumps ratio |AT"|/|AT|. Substituting 
(111 Oft into (|43|) leads after some algebra to 

K*(t w ,s) = K (t w )K 1 {s), (118) 

where s = t — t w , 

Ko(t w ) = — — 6 — , (119) 

1 - 9E{tw) 

K^s) = e~ xWs -e- x(2)s . (120) 

A factorization similar to ( I118P was found in the previous section for a general Markov 
system such that ^ 0, equation (1741) . but there it was restricted to the limit 

I AT' I I AT I or, equivalently, large values of t w . Here, it has been obtained as a 
general result for the disordered TLS in the linear response approximation. 
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The position of the maximum of K* as a function of s does not depend on t w , and 
it is given by 



1 A0> 

- A« A« 
The maximum value, for given t w , is 



where the parameter 
A« 

y = ^<l (123) 
has been introduced. In the limit t w — > (|AT| <C | AT" | ) it is 

\(2) - X(!) 

^-^ tWAW + tWAW (124) 
showing that the higher the Kovacs effect the further from a single exponential the shape 
of the relaxation. On the other hand, for t w — > oo (| AT| 3> |AT'|), 

K (t w ) <1. (125) 

Note that given the independence of s/~ from equation (I124p requires that 3> t w , 
while equation dl25[) holds for Sk <C This is due to the factorization of K*(s) and it 
seems to be also the case in other simple models [TU HE]. Between the two above limits, 
Ko(t w ) decreases monotonically, since K' Q {t w ) < for all t w , as it can be easily checked 
by direct computation. Therefore, 

^f^ = K' (t w )K l (s k )<0. (126) 

In figure [31 the Kovacs hump function K* is plotted as a function of time s for 
different values of the intermediate temperature T 1; as indicated in the insert. For all 
the curves, T = 0.43, V- = 3, V + = 7, and 7 = 0.1. The units of temperature, energy 
and time have been fixed by taking kg = 1, e = 1 and a = 1, respectively. The initial 
temperature is such that T = 0.999To. A very good agreement is found between the 
numerical solution of the master equation and the theoretical predictions, for all the 
cases considered. The curves have a maximum which is located at a fixed position, 
accurately predicted by (11211) which gives s k ~ 10.4, i.e. log 10 s k ~ 1.04. The height of 
the maximum is a decreasing function of T 1; i.e. an increasing function of the second 
temperature jump |AT'|, as implied by (11261) . 

As a consequence of the factorization property given in (11181) . 

K*(t w ,s) = K^s) 

K* max (t w ) Kfa)' 1 ) 

that is independent from t w . This is verified in figure HI where the same data as in figure 
[3] are shown, but now scaling the values of K*(t w ,s) with K^^ituf). It is seen that all 
the curves collapse as required by (11271) . 
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Figure 3. Plot of the dimensionless Kovacs function K*(t w ,s) as a function of 
the time after the second temperature jump s, for different values of the second 
temperature jump, and fixed values of the initial and final temperatures. The values 
of the parameters of the model, as well as the values of the fixed initial and final 
temperatures and the time unit used are given in the text. The symbols correspond to 
the numerical integration of the master equation while the solid lines are the theoretical 
predictions, given by (I118|) . 
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Figure 4. Plot of K*(t w , s)/K^ iax (t w ) as a function of time s, for the same values 
of the parameters as in figure [3] Time is measured in units of the TLS attempt rate 
a. Note the logarithmic scale in the vertical axis. The symbols correspond to the 
numerical integration of the master equation and the line to the analytical expression 
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Figure 5. Plot of K*(t w , s) for the same systems as in figure[3l in the long time region 
100 < s < 1000. Time s is measured in units of the TLS attempt rate a. The symbols 
correspond to the numerical integration of the master equation, while the lines are the 
theoretical result ([118JI . 
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Now, we are going to study the long time behaviour of K*(s). From (1118p - (1120p . 
it is found that for t w — > oo 

In K*(t w ,s) ~ \nK (t w ) - A (1) s , s ->■ oo , (128) 

or, taking into account the form of the long time limit of <Pe{s) (11121) . 

Ko(t w 



\nK*(t w ,s) ~ \n<j) E (s) + In 



(129) 



6W 

The difference between In K*(t w , s) and ln0 s (s) becomes independent of s for s — > oo, 
consistently with fl82|) . In figure [51 the Kovacs function is plotted, for the same values 
of the parameters as in figure [3l in the long time window 10 2 < s < 10 3 . All the curves 
are straight lines of slope — A^, as predicted by (I128p . It is also seen that for small 
values of the waiting time, corresponding in the figure to T/T\ > 1.02, the long time 
behaviour of K*(t w , s) becomes independent of t w , as it follows from (I124p . 

For high 7, namely 7 3> 25(1 + p), the parameter k defined in fl98l) behaves as 

27 

while the corresponding limits for \( 2 \ and are given by (1 114 aft and ( 1114&p . 
Moreover, these limits imply that the relaxation function </>e(t) tends to an exponential, 
and, consequently, the Kovacs effect decreases. In particular, taking into account that 
y = A (1) /A (2) -> 0, it follows from (I122p that K^ x also goes to zero. 

This is illustrated in figure [61 where the Kovacs function is plotted for two different 
high values of 7. The rest of the parameters are the same as in figure [TJ The 
temperatures values are T = 0.999T = l.OOlXi. The height of the maximum decreases 
with 7, while its width increases. The increase of the width can be understood by 
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K*(xl0 6 ) 




Figure 6. Plot of K* . for the same values of the parameters of the model of figure 
[31 for large 7, i.e. rapid barrier fluctuations. Units are the same as in the previous 
figures. The symbols correspond to the numerical integration of the master equation, 
while the lines are the theoretical expression (|118[) . Note that K* is very small. 



calculating the second derivative of K* at the maximum. A straightforward calculation 
yields 

K*"(s k ) = K {t w )K'i{s k ) = K»{t w )\W y^y {y 2 - y) , (131) 
As y — > in the limit of large 7, it follows that 

K*"(s k ) = -b^X^ 2 6 — e -- y^y , (132) 

1 - <Mm 



since is seen to go to a constant by making use of ( l91~i) . ( 111161) and ( I130p . 



Therefore, the maximum becomes flatter as 7 increases, in agreement with the behaviour 
shown in figure El 

Consider now the regime of slow barrier fluctuations, 7 <C 25(1 + p), for which k 
becomes 

«~<J(1 + /£>), (133) 

and the corresponding behaviours of \^ 2 \ and b^ are given by flll5a[) . (1115fej) . 
and (11161) . As discussed around those equations, the static disorder limit is recovered and 
the nonexponential character of the relaxation is maximal. Therefore, the Kovacs effect 
is expected to increase as 7 decreases, reaching a maximum in the static disorder limit. 
This is clearly observed in figure [TJ where K*(t w , s) is plotted for the same values of the 
other parameters as in figure El except that three small values of 7 are considered now. 
The curve corresponding to the lowest value of 7, namely 7 = 10 -4 is indistinguishable 
from the results for the static disorder. The above discussion supports the idea that the 
Kovacs effect can be understood as a measure of the non-exponential character of the 
relaxation, in the sense that the height of the peak increases and its width decreases as 
the relaxation function separates from the exponential decay. 
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Figure 7. The Kovacs function K* , for the same values of the parameters as in figure 
IH1 except that now 7 is chosen in the regime of slow fluctuations of the barrier, as 
discussed in the main text. The symbols correspond to the numerical integration of 
the master equation, while the lines are the theoretical expression (|118[) . 



Finally, the variation of the Kovacs effect with the temperature is shown in figure 
El where the dynamical disorder rate is 7 = 0.01 in all cases. The temperature jumps 
are the same as in Figs. Eland [71 T = 0.999T = l.OOITi, but several values of the final 
temperature T have been considered. The height of the maximum K^ ax as a function of 
T exhibits a non-monotonic behaviour, with a maximum for T rj 0.63. This "resonant" 
Kovacs behaviour can be easily understood. Both in the high and the low temperature 
regimes, the relaxation is almost exponential. The reason is that the rapidly fluctuating 
barrier condition 7 ^> 5(1 + p) is verified for T — > 00 and also for T — > 0, because 5 — > 
in both limits. For high temperatures, it is V_,V + <C k B T, and from (183 a\i and (183 b\i 
it follows that V- ~ v + — > a, i.e. the difference between the barriers become negligible, 
and the system is equivalent to a TLS with an average barrier. The situation in the 
low temperature regime is more subtle. In the limit ksT <C V_, V+, it is ^ v + and, 
for instance, a system of two independent TLS with those parameters should exhibit 
a strong nonexponential behaviour. But in the system with dynamical disorder, when 
both rates V- and v + tend to zero while 7 reamins constant, the fluctuations of the 
barrier are much more rapid than the transitions over it and, consequently, the system 
is again equivalent to a TLS with an average barrier. This means that K^ ax will show 
a maximum as a function of the temperature when 7 is of the order of the difference 
between the characteristic rates of the TLS with barriers V_ and V+. Then the maximum 
must occurs roughly at a the temperature T such that 

7^25(l + p). (134) 

For the values of the parameters of figure U T ^ 0.627, in very good agreement with 
the reported numerical data. 




Figure 8. The relaxation function K* for the same values of the parameters of 
the model of the previous figures, but varying the value of the final temperature T. 
The values of the ratio between the initial and intermediate temperatures have been 
kept constant, and they are given in the main text. The Kovacs effect shows a non- 
monotonic behaviour, being maximal for T w 0.63. The symbols correspond to the 
numerical integration of the master equation. For the sake of clarity, the theoretical 
lines have not been plotted. 



How good remains the linear response approximation as the T jumps are increased? 
For the range of parameters considered in this section, the agreement between the 
theoretical and the numerical integration curves remains good up to \AT\/T and 
\AT'\/T of about 0.01. Here, good means that the difference between both curves is not 
seen over the scale of the figures. If one of the jumps is further increased, multiplying 
it by a factor, the other has to be accordingly decreased, dividing it by the same factor, 
in order to keep the relative error roughly constant. 

7. Conclusions 

In this paper, the Kovacs or crossover effect has been analyzed for systems whose 
dynamics is described by a master equation. The temperature jumps have been assumed 
to be small enough so as to apply the linear response approximation. This has allowed 
to write down a general expression for the relaxation function K* describing the Kovacs 
effect, in terms of the linear relaxation function <pE arid the ratio of the temperature 
jumps. The structure of this general expression ( H3|) . is similar to the one found 
by Kovacs in his phenomenological description [5]. The main difference lies in the 
arguments of the linear relaxation function <pE- While here they are real times, in Kovacs' 
expressions they are some effective times, in the spirit of the Narayanaswami-Moynihan- 
Tool phenomenological theory of supercooled liquids j2TJ E2] E3]. Interestingly, these 
effective times can be seen to reduce to real times if the temperature jumps are 
considered to be small enough. 
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Some general properties of the function K* have been proved. They can be 
summarized as follows. The function K* is bounded, < K*(t) < 4>E{t) < 4>e{s), with 
s = t-t w . Besides, the behaviours of K*(s) in the limits |AT'| > |AT|, |AT'| < |AT|, 
and t — > oo are controlled by the behaviour of the logarithmic derivative of the direct 
relaxation function <pE{t), i.e. of the function D\{f) introduced in (jSj). If Di(t w ) has 
well defined limits for t w — > and t w — > oo, different from both zero and infinity, the 
position of the maximum Sk has also well defined limits in both cases. On the other 
hand, if -Di(O) — >■ oo, the position of the maximum Sk goes to zero for t w — > 0, i.e. 
when | AT' | ^> |AT|. Besides, if .D^oo) — > 0, Sk diverges for t w — > oo . In general, the 
amplitude of the maximum increases with the second temperature jump | AT'|. For very 
long times, K*(s) tends asymptotically to </>e(s) if the long time limit of Di vanishes. In 
this case, the direct relaxation decay, without second temperature jump, is recovered. If 
-Di(oo) ^ 0, the recovery of the direct relaxation behaviour happens in a weaker sense. 
It can only be assured that In K*(s) — ln0£(s) becomes a function of the temperature 
jumps ratio |AT'|/|AT| for large s. Then, all the curves of In K*(s) corresponding to 
different second jumps, but to the same initial and final temperatures, can be collapsed 
on the curve ln0g(s) by subtracting an adequate constant to each of them. 

The experimental observations reported by Kovacs [4, 5J are fully consistent with 
the picture discussed above for the case in which T*i(0) — > oo and .D^oo) — > 0. This is 
reasonable, since the the KWW relaxation, often found in supercooled and other slowly 
relaxing systems, is also associated to this behaviour of the logarithmic derivative of 
the relaxation funtion. An important conclusion of the analysis carried out here is 
that the Kovacs effect is still present if the temperature jumps are small, i.e. the 
experimentally observed features persist in the linear response regime. The use of the 
detailed balance condition in the analysis developed here deserves some comments. This 
condition is actually necessary only to assure that the relaxation function of the energy is 
a monotonically decreasing function of time, namely a linear combination of exponentials 
with all the coefficients being positive. But the main results of section [3] remain valid if 
detailed balance does not hold, but the right eigenvectors of the transition rates matrix 
W are still a basis for the solutions of the master equation. Then, making use of both 
the left and right eigenvectors of W, the relaxation function (pE would be again a linear 
combination of exponentials. The problem is that, for this situation, the coefficients of 
the modes are not necessarily positive. But, in those cases in which they are, <$e will 
present a monotonic decay as a function of time, as given by ( }34"1) . and all the subsequent 
results of Sees. [3] and 0] will remain valid. 

As an application of the general theory, a simple model has been considered. It 
is a two-level system with dynamical disorder. For this model, it has been possible 
to obtain an analytical closed equation for K*(t w ,s), which is found to factorize, 
K*(t w ,s) = K (t w )Ki(s). As a consequence, the position of the maximum Sk does 
not depend on t w , i.e. on the second temperature jump. This is consistent with the 
behaviour of Di, which has finite, nonzero, limits for both t — > and t — > oo. The 
independence of Sk on the second temperature jump has also been observed in other 
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models [MJ [TH]. The height of the maximum, ^ ax , of K* is an increasing function 
of | AT' | in its whole range of variation, since K' (t w ) < 0. Due to the factorization, 
all the curves corresponding to the same initial and final temperatures, but different 
second temperature jumps, collapse if they are rescaled with K^ x . The dependence 
of the Kovacs effect on the dynamical disorder has also been investigated, for fixed 
temperatures. For a rapidly fluctuating barrier, the Kovacs effect disappears, because 
-^■max ~ 0- This is logical, since our system is equivalent to a TLS with an average barrier 
in this regime, and the relaxation is exponential. On the contrary, for a slowly fluctuating 
barrier, the Kovacs effect is maximal and the static disorder limit is recovered. The 
model becomes equivalent to two independent TLS with different barriers. Finally, we 
have also analyzed the dependence of the Kovacs effect on the final temperature, keeping 
constant the relative magnitude of the jumps and the dynamical disorder. Interestingly, 
a non-monotonic dependence of ^ ax on the temperature shows up. This is a resonance 
like phenomenon, since the temperature for which K^ ax is maximal is the one at which 
two characteristic rates of the system are of the same order, as given by ( 11 34ft . 
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Appendix A. Linear response function at equilibrium 

Here the linear response function RA(t,t'), defined in (120]) . will be particularized for an 
initial equilibrium state, V e j(T). The form of the equilibrium distribution function is 
given in (l22l) . Then, it is 

V i °\t')=V e , k (T m ), (A.l) 

for if > t w . Similarly to the decomposition carried out in (|12p . the equilibrium 
distribution corresponding to T = T^> + AT can be written as 

V e (T) =V e (T {0) )+V' e . (A.2) 

It is 

WP e (T) = W {0) V' e + W'P e (T (0) ) + W'T' e = . (A.3) 

In the linear approximation being considered in this paper, the term WV' e is neglected 
in the above relation, resulting that 

W'P e (T (0) ) = -W (0) 7^ . (A.4) 

use of this relationship into ( 1201) yields 

ReAt -o = -EEE A rt%mw$v' e , k . (a. 5 ) 

i j k 
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The transition probability T t obeys the "backwards" equation [19] 

^m = j2T t {0 \mwff, (a.6) 



so that (1A.5|) is equivalent to 



i j 

d 



From and f lA~2|) it follows that 

Ki = v ^( Ti0) + AT ) ~ Ve, 3 (T {Q) ) « - - (^)i 0) ] A/3V ed (T^) , (A.8) 

where (E)^ is the average value of the energy in the equilibrium state at temperature 
T(°\ By combining (1 A. 71) and (I A. 8 1) it is obtained that 

ReAt-t') = -W^^^AiEjTMmrejiT®), (A.9) 

or, 

it^t - = (A(t)i5(f )>?> = A/?| (A(t)E(t'))i 0) , (A.10) 

where (A(t)E(t'))^ is the equilibrium time correlation function of the property A and 
the energy of the system at temperature T^°\ The general definition of these correlation 
functions is given in ( l25l) . The above equation is a particular case of the well-known 
fluctuation-dissipation theorem, relating the response to a perturbation of a system at 
equilibrium with a time correlation function. A similar derivation of (lA.lOj) can be found 
in 



Appendix B. Equilibrium time correlation functions 

Suppose a master equation whose transition rates verify the detailed balance condition, 
Wt j V eJ = W j iP eti . (B.l) 

where the equilibrium distribution has the canonical form given in (1221) . Equilibrium 
time correlation functions are defined by (see (1251) ) 

(A(t)B(t')) e = (A(t - t')B(0)) e = E E AB 3 T t ^\j)Ve, ■ (B.2) 

» i 

When detailed balance holds, the master equation can be solved by means of the 

eigenvalues and eigenvectors method [19J. Let us briefly summarize it. Consider the 
eigenproblem 

J^W^f = -X^4 a) . (B.3) 
j 
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For the sake of simplicity, it is assumed in the notation that the eigenvalue — A^ 
is non-degenerate, but the extension to the degenerate case is straightforward. All the 
eigenvalues — A^ are semi-defined negative, i.e. A^ > 0, and the only eigenvector 
corresponding to the null eigenvalue is the equilibrium distribution V e i- Then it is [19], 



V 

a ' J 



where the sum extends over all the eigenvectors of W. Substitution of this expression 
into flR2|) leads to 

(A(t)B(t')) e = J2^A l 4 a) j:B J ^e-^^ 

a i j 

= (A) e (B) e + ^'MWAfiWe-^^ . (B.5) 

as 

Here the prime to the right of the summation symbol indicates that the equilibrium 
eigenvector is excluded, and 

AA™ = A4 a) , AB ^ = B i4 a) ■ (B.6) 

i i 

Putting A = B, the expression for the time autocorrelation function is obtained, 

(A(t)A(t')) e = (A)l + (AA^) 2 e - x(a) ^ , (B.7) 

a 

and particularization for t = t' gives the second moment of the equilibrium fluctuations 
of the property A, 

(A\=(A)l + Y^(^) 2 . (B.8) 

a 

This equation shows that (AA^)) 2 is the contribution of the a-th mode to the dispersion 
(AA) 2 = (A 2 ) e — {A) 2 e . Equation fl30l) follows by particularizing (IB.7|) for A being the 
energy of the system. 
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